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Abstract This chapter reports on theoretical protocols for generating nonclassical 
states of light and mechanics. Nonclassical states are understood as squeezed states, 
entangled states or states with negative Wigner function, and the nonclassicality can 
refer either to light, to mechanics, or to both, light and mechanics. In all protocols 
nonclassicallity arises from a strong optomechanical coupling. Some protocols rely 
in addition on homodyne detection or photon counting of light. 



1 Introduction 

An outstanding goal in the field of optomechanics is to go beyond the regime of 
classical physics, and to generate nonclassical states, either in light, the mechanical 
oscillator, or involving both systems, mechanics and light. The states in which light 
and mechanical oscillators are found naturally are those with Gaussian statistics 
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with respect to measurements of position and momentum (or field quadratures in 
the case of light). The class of Gaussian states include for example thermal states of 
the mechanical mode, as well as its ground state, and on the side of light coherent 
states and vacuum. These are the sort of classical states in which optomechanical 
systems can be prepared easily. In this chapter we summarize and review means to 
go beyond this class of states, and to prepare nonclassical states of optomechanical 
systems. 

Within the family of Gaussian states those states are usually referred to as non- 
classical in which the variance of at least one of the canonical variables is reduced 
below the noise level of zero point fluctuations. In the case of a single mode, e.g. 
light or mechanics, these are squeezed states. If we are concerned with a system 
comprised of several modes, e.g. light and mechanics or two mechanical modes, 
the noise reduction can also pertain to a variance of a generalized canonical vari- 
able involving dynamical degrees of freedom of more than one mode. Squeezing of 
such a collective variable can arise in a state bearing sufficiently strong correlations 
among its constituent systems. For Gaussian states it is in fact true that this sort of 
squeezing provides a necessary and sufficient condition for the two systems to be in 
an inseparable, quantum mechanically entangled state. Nonclassicality within the 
domain of Gaussian states thus means to prepare squeezed or entangled states. 

For states exhibiting non Gaussian statistics the notion of nonclassicality is less 
clear. One generally accepted criterium is based on the Wigner phase space distri- 
bution. A state is thereby classified as non classical when its Wigner function is non 
positive. This notion of nonclassicality in fact implies for pure quantum states that 
all non Gaussian states are also non classical since every pure non Gaussian quan- 
tum state has a non positive Wigner function. For mixed states the same is not true. 
Under realistic conditions the state of optomechanical systems will necessarily be a 
statistical mixture such that the preparation and verification of states with a non pos- 
itive Wigner distribution poses a formidable challenge. Paradigmatic states of this 
kind will be states which are close to eigenenergy (Fock) states of the mechanical 
system. 

Optomechanical systems present a promising and versatile platform for creation 
and verification of either sort of nonclassical states. Squeezed and entangled Gaus- 
sian states are in principle achievable with the strong, linearized form of the radia- 
tion pressure interaction, or might be conditionally prepared and verified by means 
of homodyne detection of light. These are all "Gaussian tools" which conserve the 
Gaussian character of the overall state, but are sufficient to steer the system towards 
Gaussian non classical states. In order to prepare non Gaussian states, possibly with 
negative Wigner function, the toolbox has to be enlarged in order to encompass 
also some non Gaussian instrument. This can be achieved either by driving the op- 
tomechanical system with a non Gaussian state of light, such as a single photon 
state, or by preparing states conditioned on a photon counting event. Ultimately 
the radiation pressure interaction itself is a nonlinear interaction (cubic in annihi- 
lation/creation operators) and therefore does in principle generate non Gausssian 
states for sufficiently strong coupling go at the single photon level. Quite generally 
one can state that some sort of strong coupling condition has to be fulfilled in any 
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protocol for achieving a nonclassical state. Fulfilling the respective strong coupling 
condition is thus the experimental challenge on the route towards nonclassicality in 
optomechanics. 

In the following we will present a selection of strategies aiming at the preparation 
of nonclassical states. In Sec.[2]we review ideas of using an optomechanical cavity 
as source of squeezed and entangled light. Central to this approach is the fact that 
the radiation pressure provides an effective Kerr nonlinearity for the cavity, which is 
well known to be able to generate squeezing of light. In Sec.[3]we discuss nonclas- 
sical states of the mechanical mode. This involves e.g. the preparation of squeezed 
states as well as non Gaussian states via state transfer form light, continuous mea- 
surement in a nonlinearly coupled optomechanical system, or interaction with single 
photons and photon counting. Sec.|4]is devoted to nonclassical states involving both 
systems, light and mechanics, and summarizes ideas to prepare the optomechanical 
system in an entangled states, either in steady state under continuous wave driving 
fields, or via interaction with pulsed light. 



2 Non-classical states of light 
2.1 Ponderomotive squeezing 

One of the first predictions of quantum effects in cavity optomechanical system con- 
cerned ponderomotive squeezing [1 2], i.e., the possibility to generate quadrature- 
squeezed light at the cavity output due to the radiation pressure interaction of the 
cavity mode with a vibrating resonator. The mechanical element is shifted propor- 
tionally to the intracavity intensity, and consequently the optical path inside the 
cavity depends upon such intensity. Therefore the optomechanical system behaves 
similarly to a cavity filled with a nonlinear Kerr medium. This can be seen also by 
inserting the formal solution of the time evolution of the mechanical displacement 
x(t) into the Quantum Langevin equation (QLE) for the cavity field annihilation 
operator d(t), 

-+/a> opt (0) d + / dsXAf(t—s) [ihG 2 d(t)d'(s)d(s) + iGd(t)^(s)] +y/Ka in {t) 



where a,„(f) is the driving field (including the vacuum field) and 



(1) 



d(0 e~ im e -rM'/2 

%M{t)= I IT- 7— 2 1 — - — r = ~- sm£l M t (2) 

2% m ef{ (Ofa -CO 2 - ir M CO) m eff i2 M 



is the mechanical susceptibility (here Qm = y Q^—T^/4). Eq. (li shows that 
the optomechanical coupling acts as a Kerr nonlinearity on the cavity field, but 
with two important differences: i) the effective nonlinearity is delayed by a time 
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depending upon the dynamics of the mechanical element; ii) the optomechanical 
interaction transmits mechanical thermal noise %(t) to the cavity field, causing 
fluctuations of its frequency. When the mechanical oscillator is fast enough, i.e., 
we look at low frequencies CO <C Qm, the mechanical response is instantaneous, 
%u{t) — 5(f)/m e ff^M' an< ^ tn e nonlinear term becomes indistinguishable from a 
Kerr term, with an effective nonlinear coefficient = nG 2 /m^fQ^. 

It is known that when a cavity containing a Kerr medium is driven by an intense 
laser, one gets appreciable squeezing in the spectrum of quadrature fluctuations at 
the cavity output [3 1. The above analogy therefore suggests that a strongly driven op- 
tomechanical cavity will also be able to produce quadrature squeezing at its output, 
provided that optomechanical coupling predominates over the detrimental effect of 
thermal noise 0] |2] . 

We show this fact by starting from the Fourier-transformed linearized QLE for 
the fluctuations around the classical steady state 

m eff (Q^- co 2 - icor M ) x(co) = hGa s 8X(co) + t;(co), (3) 
(|-"») SX(co) = -A8Y{co) + ^8X in (co), (4) 
(| -«'<») SY{(o) = A8X(co) + Ga s x(co) + ^8Y in (co), (5) 

where we have chosen the phase reference so that the stationary amplitude of the 
intracavity field a s is real, SX = 8a + Sa^ [8Y = —i (5a — Sa')] is the amplitude 
(phase) quadrature of the field fluctuations, and 8X' n and 8Y'" are the correspond- 
ing quadratures of the vacuum input field. The output quadrature noise spectra are 
obtained solving Eqs. ([3])-(j5]l, and by using input-output relations Q, the vacuum 
input noise spectra S^(co) = Sy (©) = 1, and the fluctuation-dissipation theorem for 
the thermal spectrum (co) = ftcaZMfWeffCoth (Ticq/UcbT). 

The output light is squeezed at phase (j) when the corresponding noise spec- 
trum is below the shot-noise limit, S™'(cd) < 1, where S"^"(co) = Sf" (to) cos 2 (j) + 
S Y "'(co) sin 2 +SJy(a))sin20, and the amplitude and phase noise spectra S'^"(co) 
and S" l "(co) satisfy the Heisenberg uncertainty theorem S™' (co)S'y" (co) > 1 + 

[Sjl'y (o)] 2 [4]. However, rather than looking at the noise spectrum at a fixed phase of 
the field, one usually performs an optimization and considers, for every frequency 
CO, the field phase <p opt (co) possessing the minimum noise spectrum, defining in this 
way the optimal squeezing spectrum, 

S opt (co) = mmST(co) 

2Sr(co)sr(co)-2[ST Y (co)} 2 ... 

= 1 (6) 

S° x ut (co) + S°y u, {co) + \j[Sf{G>) - S Y ut (co)} 2 + 4 [S'^(co)} 2 
The frequency-dependent optimal phase is correspondingly given by 
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<l>opt{(o) = -arctan 



2S$(a>) 



ST(co)- 



-S'f(co) 



(7) 



Since there is no benefit for choosing a nonzero detuning, we restrict to the resonant 
case 4=0, which is always stable and where expressions are simpler. One gets 



rout 



(co) 



1. 



STr(co) 



KhG 2 ajRz{xM{co)} 



K 2 /4 + C0 2 



(8) 
(9) 



where 



Sr(C0) 



KhG 2 a 2 lm{x M (co)} 



74- 



KhG 2 a 2 \m{% M {<o)} , / hto \ 



K 2 /4- 



coth 



\2k B T 



(10) 



Inserting Eqs. ([8]l-(|9]) into Eq. |6]l one sees that the strongest squeezing is obtained 
when the two limits S r (co) <C 1 and [S^y(o))] 2 3> 1 are simultaneously satisfied. 
These conditions are already suggested by Eq. (JTJ): S r ((o) <C 1 means that thermal 
noise is negligible, which occurs at low temperatures and small mechanical damping 
I m {/frw(ft>)}> i-e-, large mechanical quality factor Q; [5^y(a))] 2 3> 1 means large 
radiation pressure, achieved at large intracavity field and small mass. Ponderomotive 
squeezing is therefore attained when 



Sr(C0) 



P m co L ,? 2 Q 
m eff c 2 ^ n t h 



»1, 



(11) 



and in this ideal limit S opt (co) ~ [5^(©)] ->• 0, and (p opt ((o) ~ - arctan [2/5^(ffl)] /2 - 
0. Since the field quadrature 8X" ut at = is just at the shot-noise limit (see 
Eq. ([H}), one has that squeezing is achieved only within a narrow interval for the ho- 
modyne phase around <j) opt (to), of width ~ 2 1 § opt (co) | ~ arctan 1 2 /S'^f ( 03 ) I ■ This ex- 
treme phase dependence is a general and well-known property of quantum squeez- 
ing, which is due to the Heisenberg principle: the width of the interval of quadrature 
phases with noise below the shot-noise limit is inversely proportional to the amount 
of achievable squeezing. 

Sopt(co) and the corresponding optimal phase § opt (co) at which best squeezing 
is attained for each co, are plotted in Fig. [TJfor a realistic set of parameter values 



(see figure caption). S opt (co) 



is below the shot-noise limit whenever S£y((o) 



^0 



(see Eqs. (|6])-(|9]l), and one gets significant squeezing at low frequencies, well below 
the mechanical resonance, where the optomechanical cavity becomes fully equiv- 
alent to a Kerr medium, as witnessed also by the fact that (p op ,(co) is constant in 
this frequency band. This equivalence is lost close to and above the mechanical 
resonance, where squeezing vanishes because Re{^(o)} ~ ^xyi 03 ) ~ 0, and the 
optimal phase shows a large variation. 

The present treatment neglects technical limitations: in particular it assumes the 
ideal situation of a one-sided cavity, where there is no cavity loss because all photons 
transmitted by the input-output mirror are collected by the output mode. We have 
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Frequency ra(Hz) Frequency to(Hz) 

Fig. 1 Optimal spectrum of squeezing in dB S„ pl (a), and the corresponding optimal quadrature 
phase f opt (b), versus frequency in the case of a cavity with bandwidth K = 1 MHz, length L = 1 
cm, driven by a laser at 1064 ran and with input power 3?i„ = 10 mW. The mechanical resonator 
has Qyi/lK = 1 MHz, mass m e g =100 ng, quality factor Q = 10 4 , and temperature I=4K. 



also ignored laser phase noise which is typically non-negligible at low frequencies 
where ponderomotive squeezing is significant. In current experimental schemes both 
cavity losses and laser phase noise play a relevant role and in fact, ponderomotive 
squeezing has not been experimentally demonstrated yet. However, the results above 
show that when these technical limitations will be solved, cavity optomechanical 
systems may become a valid alternative to traditional sources of squeezing such as 
parametric amplifiers and Kerr media. 



2.2 EPR correlated beams of light 

Optomechanical cavities provide a source not only of squeezed light but also of 
entangled light, as we will explain in the following. By means of spectral filters, 
the continuous wave field emerging from the cavity can be split in many travel- 
ing modes thus offering the option of producing and manipulating a multipartite 
system Q. In particular we focus on detecting the first two motional sidebands 
at frequencies © op t ± and show that they posses quantum correlations of the 
Einstein-Podolsky-Rosen type |[6l . 

Using the well-known input-output fields connection d"'"(t) — \/2Ka{t) — a" l (t), 
the output mode can be split in A 7 * independent optical modes by frequency selection 
with a proper choice of a causal filter function: 

&T (t) = f__ dsg k (t - s)a out (s), k=l,...N,, (12) 

where gk(s) is the causal filter function defining the k-th output mode. The annihila- 



te, 



tion operators describe N independent optical modes when 
which is fulfilled when /J° dsgj^ygkis) — 8m, i.e., the N filter functions gk(t) form 



if{t),a° k ut (tf 
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Fig. 2 a) Logarithmic negativity of Stoke-Antistokes output modes when Q\ = —Qm while Qj is 
varied around Qm- The inverse bandwidth is kept constant at e = I07t. b) Temperature robustness 
of bipartite entanglement of output modes at ±Qm computed for short (e = 10?r, dashed line) 
and long (e = IOOtt, solid line) detection times, c) The bipartite Stoke-Antistoke entanglement 
shows improvement and eventually saturates with increasing the integration time. Parameters are 
Qyi/Xn = 10 MHz, Q = 10 s , mass m e ff = 50 ng, cavity of length L = 1 mm with finesse & = 
2 x 10 4 , detuning A = Qm, input power P m = 30 mW at 810 nm, and temperature T = 0.4 K, 
yielding go = 0.43 kHz, g = 0.4112m, a cavity bandwidth K = 0.7512m, and a thermal occupation 
of n t h — 833. 



an orthonormal set of square-integrable functions in [0,°°). As an example of a set 
of functions that qualify as casual filters we take 

8(r)-e(r-T) _ inkt 

gk{t) = 7= e k , (13) 

v T 

(0 denotes the Heavy side step function) provided that and T satisfy the condition 
£2 j ■ — £2k — ^-p for integer p. Such filtering is seen as a simple frequency integration 
around £2^ of bandwidth ~ \ Jx (the inverse of the time integration window). 

For characterization of entanglement one can compute the stationary (2A^ + 2) x 
(2/V + 2) correlation matrix of the output modes defined as 

V" u, (t) = - Uf(t)u'f(t) + u'' ut (tX ut (t)) , (14) 

where 

U ""'( t ) = (q(t) 7 p( t ),xr( t ),Yr(t), . . . fsrwfwf 

is the vector formed by the mechanical position and momentum fluctuations and by 
the amplitude (t£ ut (t) = [af"(f ) + af{t ) f ]/V2), and phase (tg ut (t) = [a° k ut (f) - 
d% ut (t) 7 } / iy/2) quadratures of the output modes. 

We are now in position to analyze the quantum correlations between two output 
modes with the same bandwidth T _1 and central frequencies £2\ and Cli- As a mea- 
sure for entanglement we apply the logarithmic negativity E jy to the covariance ma- 
trix of the two optical modes. It is defined as E , ,y = max[0, — ln2rj~], where r\~ = 
2 -i/2 



1 /2 

E(V) - [Z{V f -4detV] 1/2 , with E(V) = detV m + detV c - 2detV m 
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and we have used the 2x2 block form of the covariance matrix 

Therefore, a Gaussian state is entangled if and only if 77 < 1/2, which is equiva- 
lent to Simon's necessary and sufficient entanglement non-positive partial transpose 
criterion for Gaussian states, which can be written as 4det V < E — 1/4. 

The resulting quantum correlations among the upper and the lower sideband in 
the continuous wave output field are illustrated in Fig. [2] We plot the interesting 
and not unexpected behavior of E jf as a function of central detection frequency Q 
in Fig j2^, with the mirror reservoir temperature in Fig|2j3 and with the scaled time 
integration window e = X2mT in Figj2j;. The conclusion of Fig|2^ is that indeed scat- 
tering off the mirror can produce good Stokes-Antistokes entanglement which can 
be optimized at the cavity output by properly adjusting the detection window. More- 
over, further optimization is possible via an integration time increase as suggested 
by Fig|2j;. The temperature behavior plotted in Fig|2j5 shows very good robustness 
of the mirror-scattered entangled beams that suggests this mechanism of producing 
EPR entangled photons as a possible alternative to parametric oscillators. 



3 Non-classical states of mechanics 

3.1 State transfer 

For a massive macroscopic mechanical resonator, just as in the case of a light field, 
the signature of quantum can be indicated in a first step by the ability of engineering 
of a squeezed state. Such a state would also be useful in ultrahigh precision mea- 
surements or detection of gravitational waves and has been experimentally proven 
in only one instance for a nonlinear Duffing resonator [7]. Numerous proposals exist 
and can be categorized as i) direct: modulated drive in optomechanical settings with 
or without feedback loop 1 8 9l [T0l[TT1l . and ii) indirect: mapping a squeezed state of 
light or atoms onto the resonator, coupling to a cavity with atomic medium within 
fl2l . coupling to a Cooper pair box [ 13 1 or a superconducting quantum interference 
loop lfT4l IT31 . In the following we take the example of state transfer in a pure op- 
tomechanical setup where laser cooling of a mirror/membrane via a strong laser is 
accompanied by squeezing transfer from a squeezed vacuum second input light field 
fl6l . While the concept is straightforward it is of interest to answer a few practical 
questions such as: i) what is the resonance condition for optimal squeezing transfer 
and how does a frequency mismatch affect the squeezing transfer efficiency, ii) what 
is the optimal transfer, iii) how large should the cavity finesse should be for optimal 
transfer etc. 

To this purpose we assume the typical membrane in the middle optomechanical 
system depicted in Fig. [3^, where the laser Antistokes sideband, the cavity and the 
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Fig. 3 a) Logarithmic negativity of Stoke-Antistokes output modes when Q\ = —£2m while Q.^ is 
varied around Q.yy. The inverse bandwidth is kept constant at e = lOn. b) Temperature robustness 
of bipartite entanglement of output modes at ±Qm computed for short (e = 107T, dashed line) and 
long (e = IOOtt, solid line) detection times, c) The bipartite Stoke-Antistoke entanglement shows 
improvement and eventually saturates with increasing the integration time. 



squeezing central frequency are resonant (shown in [3}}). The input squeezed light 
operators have the following correlations 

(c m (t + T)c ta (0) = f ^2 (bye-W +b x e- b ^) (16) 

<4(* = f ^% (bye-^-he-^) . (17) 

y x 

The noise operators are written in a frame rotating at 0) s and satisfy the canon- 
ical commutation relation [c- m (t),c^ n (t')] — 5{t — t'). Parameters Af and M deter- 
mine the degree of squeezing, while b x and b y define the squeezing bandwidth. 
For pure squeezing there are only two independent parameters, as in this case 
\M\ 2 = N(N+ 1) and b y = b xy /2 (N + |M|) + 1. 

Following the standard linearized quantum Langevin equations approach for op- 
tomechanics, we first identify two conditions for optimal squeezing: i) A = £2m, 
meaning that we require continuous laser cooling in the resolved sideband regime 
and ii) A s = —£2m so that the squeezing spectrum is centered around the cavity 
frequency. Then we look at the variances of the generalized quadrature operator 

dl 9 {t) = ^{e'n{t)+e-'n\t)) ) (18) 

which for <p = is the usual position operator q(t) and for <p = —n/2 is the mo- 
mentum operator /3(f), both taking in a rotating frame at frequency 12m- In the limit 
of squeezed white noise the quadrature correlations take a simple form 
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(SX^SX^t)) = (N+ l --Re{Me 21 ^ +g (%+ *) ■ (19) 

The first term in the right hand side comes from the squeezing properties of the 
squeezed input vacuum while the second term is the residual occupancy after laser 
cooling. In view of this equation a successful squeezed mechanical state prepara- 
tion automatically requires close to ground state cooling. One can follow this in 
Fig(3]: where cooling close to ground state of an initially thermal mechanical state 
is performed by the cooling laser, and subsequently squeezing of a quadrature is 
achieved via the squeezed vacuum. In Fig|3jl we suggest that although ground state 
cooling is not easy to achieve, one can at least in principle (for example by complete 
tomography of the mechanical state) see a quadrature squashing effect. 

To answer a practical question, when the squeezing is not white, fulfilling the 
resonance condition A s = —£2m is important. The deviations allowed for the fre- 
quency mismatch are smaller than the width of the cooling sideband, i.e. r e g. A 
second question is the effect of the finite width of squeezing. In general there is an 
optimal squeezing bandwidth for which the transfer from light to membrane is max- 
imized, but in the resolved sideband limit where Qm 3> K - the finite bandwidth result 
does not differ much from the infinite bandwidth limit result. For a large bandwidth 
which fully covers the motional sidebands, b x ^> Q.yi, the membrane sees only white 
squeezed input noise, whereas for smaller bandwidth, the crucial question is whether 
the squeezed input will touch the heating sideband or not. For a high-finesse cavity, 
the width is not a big issue, since the heating sideband is anyway weak, whereas for 
a bad cavity the squeezing transfer is much improved for an optimal, finite band- 
width where the strong heating sideband is avoided. 



3.2 Continuous measurements of mechanical oscillators. 

Another method for creating nonclassical states exploits the possibility to condition- 
ally prepare states of the mechanical oscillator via measuring the output field of the 
optomechanical system. The coupling of a mechanical resonator to an optical cavity 
field enables an indirect continuous monitoring of the mechanical motion by a direct 
phase dependent measurement of the field leaving the cavity. The standard radiation 
pressure coupling is linear in the displacement of the mechanical resonator and thus 
enables a continuous measurement of displacement. If we wish to monitor the en- 
ergy (or phonon number) of the mechanical resonator, however, we need to find 
an interaction Hamiltonian that is quadratic in the displacement. Such interactions 
can occur in a number of ways ifTTl [18). We begin with the case of displacement 
measurements. 

The standard linearised opto-mechanical coupling Hamiltonian 

H = tiAa^a + h£2 M b f b -hg{a + a^){b + b r ) (20) 
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As the interaction part of this Hamiltonian commutes with the (dimensionless) me- 
chanical displacement operator, q = i (b + V), in principle this model can be con- 
figured as a measurement of the displacement. However as q does not commute 
with the free mechanical Hamiltonian, this is not a strict QND measurement lfT9ll . 
Nonetheless, for a rapidly damped cavity, we can effect can approximate QND read- 
out of the mechanical displacement provided the coupling constant g can be turned 
on and off sufficiently fast. This can be achieved by using a pulsed coherent driving 
field on the cavity ll20ll . 

If we include the damping of both the cavity and the mechanics, we obtain the 
quantum stochastic differential equations, 

da K ~ /— 

— = —iAa — —a + igX + y fca,„ (21) 
at 2 

— = -i£2 M b - r ^ L b + ig(a + a T ) + \Zf^b in (22) 
at 2 

where we assume that the input noise to the cavity is vacuum, so that the only 
non zero correlation function for the cavity noise is (di n (t)dj n (f)) = 8(t — t') but 
that the input noise to the mechanical resonator is thermal, {bi n {t)b] n {t')) = («th + 
l)5(f — ?'). We expect that a good measurement will occur when the cavity field is 
rapidly damped so that it is slaved to the mechanical degree of freedom. We can 
then adiabatically eliminate the cavity degree of freedom by setting to zero the right 
hand side of Eq.pT) and formally solving for the operator a, 



y/A 2 + K 2 /A iA + K/2 

where tan =2A/k. The actual output field from the cavity is related to the field 
inside by the input/output relation, a uut = y/xa — &i n , so that 

a m {t) - ^5 '[ m+e- 2i ^ in (t) (24) 

Clearly this indicates that we need to measure a particular quadrature of the output 
field (for example by homodyne or heterodyne detection) and that the added noise is 
vacuum noise. The optimal transfer is obtained on resonance. A fast, impulsive read- 
out of the mechanical resonator's displacement may be made by injecting a coherent 
pulse into the cavity and subjecting the output pulse to a homodyne measurement 

Eo). 

If we wish to measure the energy of a mechanical resonator we must find an inter- 
action hamiltonian that is at least quadratic in the mechanical amplitude. A number 
of schemes have been proposed, including trapped atoms in a standing wave ETl 
and a nanomechanical resonator coupled to a single superconducting junction dipole 
in the dispersive regime [22]. In opto-mechanics a dielectric membrane placed at the 
antinode of a cavity standing wave shifts the cavity frequency proportional to the 
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square of the mechanical displacement of the membrane from equilibrium ifTTl . A 
similar interaction arises for an optically levitated particle in a standing wave fllfl . 
The interaction Hamiltonian in this case takes the form 

H = h(O op ,a i a + hn M Pb + h(e*ae imL ' + e c <zV ,£B ' f ) + -G 2 d^a(h + b'f, (25) 
where 



h d 2 CO c {x) 

<J2 



2vm 



(26) 



x=x 

and where we have included a coherent driving field with amplitude £. As usual 
we expand the interaction around the steady state cavity field. After the rotating 
wave approximation, the effective Hamiltonian in the interaction picture may then 
be written as 

H I =^ X (S + & f )b% (27) 

where % — IGiO®. 



The interaction in Eq.(27i describes a displacement of the cavity field propor- 
tional to the number of vibrational quanta in the mechanical resonator. The average 
steady-state displacement is given by Jn/fc, where h is the mean phonon number 
operator for the mechanical oscillator, K is the cavity line-width. If we continu- 
ously monitor the output field amplitude from the cavity via homodyne detection 
this scheme can in principle enable a continuous monitoring of the mechanical vi- 
brational energy, and phonon number jumps [23 1. 

Under continuous homodyne measurement of this quadrature, the system is gov- 
erned by the following stochastic master equation (SME): 

dp = - - [Hi , p ] dt + y(n, h + 1 ) 9 [b] p dt + yn th 9 [P ] p dt + k3> [a] p dt + \fKdW Jt?[ae~ 

where 2$[c]p = cpc t — c' cp /2 — pc^c/2 and J^f[c]p = cp+pc^— Tr(cp + pc^) is 
the measurement super-operator, y and K are the respective mechanical and cavity 
damping rates. 

In figure H] we show a numerical integration of the stochastic master equation 
with K = NryN for three cases: X 2 / K = % 2 /k = 10yn,/„ and % 2 / K = 10 2 Yh th . 
We start with the mechanics in the ground state, and the bath temperature is set at 
n,i, = 0.5. The first case, Fig. |4^), does not satisfy the fast-measurement condition 
and therefore does not resolve quantum jumps in the phonon number. The second 
case, Fig. (J4j?), is on the border of the fast-measurement regime for n ~ 1 and shows 
some jump-like behaviour in the phonon number. The third case, Fig. (|4j;), strongly 
satisfies the fast-measurement condition for low phonon numbers and shows well- 
resolved quantum jumps in spite of being deeply within the weak coupling regime 
with x/K = 1CT 1 . 

For jump-like behaviour to arise in the weak coupling limit, the adiabatic con- 
dition is not sufficient. In this regime, analysis shows that the rate of information 
acquisition about the phonon number is proportional to x 2 / K. As in the strong cou- 
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Fig. 4 The evolution of the conditional average phonon number with parameters K = 10 4 yn f /, and: 
(a) X 2 l K = 7«rt> (b) X 2 l K = 10y/V, (c) £ 2 / K = 10 2 y« f /,. Jump-like behaviour occurs only when 
X 2 / K ^> y\n t f,{n + 1) + (fiii, + !)»], where n is the phonon number. 



pling case, this measurement rate must dominate the thermalisation rate in order for 
quantum jumps to arise. Thus, in addition to being in the adiabatic limit, the weak 
coupling regime requires % 2 / K 3> jN. 

The state of the mechanical system conditioned on the measurement of light 
will be in a Non-Gaussian state. This is due to the non-linear interaction introduced 
in ( |25| l. Another way of achieving a Non-Gaussian state exploits the nonlinearity 
provided in photon counting, as will be detailed in the next section. 



3.3 Non-Gaussian state via interaction with single photons and 
photon counting 

In 1935 Schrodinger pointed out that according to quantum mechanics even macro- 
scopic systems can be in superposition states |24|. The interference effects, charac- 
teristic of quantum mechanics, are expected to be hard to detect due to environment 
induced decoherence (25] . Nevertheless there have been several proposals on how 
to create and observe macroscopic superpositions in various physical systems. See 
references ll26l l27l l28l for some of the first proposals. There have also been experi- 
ments on superposition states in superconducting and piezoelectric devices ll29l[30l 
and on interference with fullerene ff3~TTl and other large molecules. One long-term 
motivation for this kind of experiment is the question of whether unconventional 
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decoherence processes such as gravitationally induced decoherence or spontaneous 
wave-function collapse [32, 33, 34, 35, 36 1 take place. 

In this section a scheme is analyzed that is close in spirit to Schrodinger's origi- 
nal discussion. A small quantum system (a photon) is coupled to a large system (a 
mirror) such that they become entangled 11371 . The system consists of a Michelson 
interferometer in which one arm has a tiny moveable mirror. The radiation pressure 
of a single photon is used to displace the tiny mirror. The initial superposition of 
the photon being in either arm causes the system to evolve into a superposition of 
states corresponding to two distinct spatial locations of the mirror. A high-finesse 
cavity is used to enhance the interaction between the single photon and the mirror. 
The interference of the photon upon exiting the interferometer allows one to study 
the creation of coherent superposition states periodic with the motion of the mirror. 

Consider the setup shown in Fig. [5] consisting of a Michelson interferometer that 
has a cavity in each arm. In the cavity in arm A one of the mirrors is very small 
and attached to a micromechanical oscillator. While the photon is inside the cavity, 
it exerts a radiation pressure force on the small mirror. We will be interested in the 
regime where the period of the mirror's motion is much longer than the roundtrip 
time of the photon inside the cavity, and where the amplitude of the mirror's motion 
is very small compared to the cavity length. Under these conditions, the system can 
be described by the standard optomechanical Hamiltonian 



H = (O opt a J <a + n M Pb-g a' ! <a(b + b' f ). (29) 

To start with, let us suppose that initially the photon is in a superposition of being 
in either arm A or B, and the mirror is in a coherent state 18) = gH^I 2 / 2 y Jt\ n \ 
where \n) are the eigenstates of the harmonic oscillator. Then the initial state is 

|^(0)) = ^(|0) A |l) fi + |l) A |0) B )|j8). (30) 

After a time t the state of the system will be given by fiOlHTI 

Ml)) = ^=e-^'{\0) A \l) B \l5e- n ^) 

+ e inHn M t-smQ Mt )^ A ^ B ^ e -i£i Mt + r] ^ _ e -i£2 M t^^ (31) 

where 77 = go /Qm- In the second term on the right hand side the motion of the 
mirror is altered by the radiation pressure of the photon in cavity A. The parame- 
ter 7] quantifies the displacement of the mirror in units of the size of the coherent 
state wavepacket. In the presence of the photon the mirror oscillates around a new 
equilibrium position determined by the driving force. 

The maximum possible interference visibility for the photon is given by twice 
the modulus of the off -diagonal element of the photon's reduced density matrix. By 



tracing over the mirror one finds from Eq. ( 3 1 1 that the off -diagonal element has the 
form 
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Fig. 5 Scheme: a Michelson interferometer for a single photon with in each arm a high-finesse 
cavity. The cavity in arm A has a very small end mirror mounted on a micro-mechanical oscillator. 
The single photon enters the interferometer via a polarizing beam splitter (PBS) followed by a X/4 
wave plate. This is an optical trick to allow detection of the photon leaking out of the interferometer 
at a later time on detector Dl or D2. If the input photon is considered to be in arm A, the motion of 
the small mirror is effected by its radiation pressure. If the input photon is considered to be in arm 
B, the motion of the mirror is undisturbed. The interferometer, based on the 50%/50% beamsplitter 
(BS), leads to the entanglement between the photon being in arm A or in arm B and the state of the 
mirror. 



where Im denotes the imaginary part. The first factor is the modulus, obtaining its 
minimum value after half a period at f = tz/Qm, when the mirror is at its maximum 
displacement. The second factor gives the phase, which is identical to that obtained 
classically due to the varying length of the cavity. 



For general t the phase in Eq. (32i depends on j3, i.e. the initial condition of 
the mirror. However, the effect of the initial condition averages out after every full 
period. 

In the absence of decoherence, after a full period, t = 2k/Qm, the system is in 
the state ( |0)a | l)s + e rr?_27r |l)^|0)^)|j3), such that the mirror is again completely 
disentangled from the photon. Full interference can be observed if the photon is de- 
tected at that moment. If the environment of the mirror "remembers" that the mirror 
has moved, then, even after a full period, the photon will still be entangled with 
the mirror's environment, and thus the interference for the photon will be reduced. 
Therefore the setup can be used to measure the decoherence of the mirror. 

In practice the mirror attached to a mechanical-resonator will be in a thermal 
state, which can be written as a mixture of coherent states |j3) with a Gaussian 
probability distribution (l/7rn t h)e _ '^' 2,/ " lh , where n t h is the mean thermal number of 
excitations, n t h = 1/ [e tlQ,M l kT — 1). If one wants to determine the expected interfer- 
ence visibility of the photon at a time t for an initial mirror state which is thermal, 



one therefore has to average the off-diagonal element Eq. ( 32 1 over j3 with this dis- 
tribution. The result is 
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_ e -Tl 2 {2n th + l){l-cosn M t) e iri 2 (n M t-smn M t) _ 



(33) 



As a consequence of the averaging of the f5 -dependent phase in Eq. (32i, the off- 
diagonal element now decays on a timescale l/(TJ^M-v/"th) a ^ el t — Q, i.e. very 
fast for the realistic case of large n t h. However, remarkably it still exhibits a revival 
l27l at t = 27z/Qm, when photon and mirror become disentangled and the phase 
is independent of j3, such that the phase averaging does not reduce the visibility. 
Figure [6] shows the time evolution of the visibility of the photon over one period of 
the mirror's motion for T} = 1 and temperatures T of 1 mK, 100 jJ.K and 10 /iK. 

The magnitude of the revival is reduced by any decoherence of the mirror. Fur- 
thermore the revival will also be reduced due to nonlinear terms in the mechani- 
cal oscillator. However since we will only consider extremely small displacements 
around the equilibrium position we assume that nonlinear effects can be ignored. 

The revival demonstrates the coherence of the superposition state that exists at 
intermediate times. For r\ 2 > 1 the state of the system is a superposition of two 
distinct positions of the mirror. More precisely, for a thermal mirror state, the state 
of the system is a mixture of such superpositions. However, this affects neither the 
fundamentally non-classical character of the state nor, as we have seen, the existence 
of the revival after a full period. We now discuss the experimental requirements for 
achieving such a superposition state and observing its recoherence at t = 2k/£2m- 




Fig. 6 Time evolution of the interference visibility of the photon over one period of the mirror's 
motion for t] = 1 and temperatures T of 1 mK, 100 /xK and 10 /iK. The visibility decays very fast 
after t = 0, but in the absence of decoherence there is a sharp revival of the visibility after a full 
period (2k). The width of each peak scales like 1 / yT. 



Firstly, it is required that r\ 2 > 1, which physically means that the momentum 
kick imparted by the photon to the mirror has to be larger than the initial quantum 
momentum uncertainty of the mirror. Let N denote the number of roundtrips of the 
photon in the cavity during one period of the mirror's motion, such that 2NL/c = 
2k/£2m- This allows us to rewrite the condition rj 2 > 1 as 

2M3L >i, (34) 
7Zcm e ffX 2 ~ 
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where A is the wavelength of the light. The factors entering Eq. (34i are not all 
independent. The achievable N, which is determined by the quality of the mirrors, 
and the minimum possible mirror size depend strongly on A. The mirror's lateral 
dimensions should be an order of magnitude larger than A to limit diffraction and 
avoid geometrical losses. The minimum possible thickness of the mirror generally 
depends on the wavelength as well in order to achieve sufficiently low transmission. 



Eq. (34 1 allows one to compare the viability of different wavelength ranges. 
While the highest values for N are achievable for microwaves (up to 10 10 ), this 
is counteracted by their long wavelengths (of order cm). On the other hand there are 
no good mirrors for highly energetic photons. The optical regime seems optimal. In 
the following estimates we will consider a A around 630 nm. 

The cavity mode needs to have a very narrow focus on the tiny mirror, which 
requires the other cavity end mirror to be large due to beam divergence. The maxi- 
mum cavity length is therefore limited by the difficulty of making large high quality 
curved mirrors. In fact from simulations it follows that indeed the surface quality of 
the large curved mirror is likely to be the most challenging component of the setup 
Il42ll43l . Here we consider a cavity length of 5 cm, and a small mirror size of 10 x 
10x10 microns, leading to a mass of order 5 x 10~ 12 kg. 

One possible path for the fabrication of such a small mirror on a good mechanical 
oscillator is to coat a silicon cantilever with alternating layers of Si02 and a metal 
oxide such as Ta20s by sputtering deposition. The best current mirrors are made 
in this way. Recently such mirrors have also been produce on Silicon Nitride cross 
resonators which have excellent mechanical properties [44 1 



For the above dimensions the condition ( 34 1 is satisfied for N = 5.6 x 10 6 . There- 
fore a photon loss per reflection not larger than 3 x 10~ 7 is needed, which is about a 
factor of 4 below the best reported values for such mirrors [45], and for a transmis- 
sion of 10~ 7 , which is consistent with the quoted mirror thickness [46 1 . For these 
values, about 1 % of the photons are still left in the cavity after a full period of the 
mirror. Coupling into a high-finesse cavity with a tight focus will require carefully 
designed incoupling optics. For the above values of N and L one obtains a frequency 
£2m = 2n x 500 Hz. This leads to a quantum uncertainty of order 10~ 13 m, which 
for r\ 2 ~ 1 corresponds to the displacement in the superposition. 

Secondly, the requirement of observing the revival after a full period of the mir- 
ror's motion puts a bound on the acceptable environmental decoherence. To estimate 
the expected decoherence we model the mirror's environment by an (Ohmic) bath of 
harmonic oscillators. Applying the analysis of 11471 one then finds that off-diagonal 
elements between different mirror positions decay with a factor 



exp 



ry\kTm e fi(Ax) 2 sin 12m? cosClut \ 

X2m / 



h 2 



(35) 



where I~m is the rate of energy dissipation for the mechanical oscillator, T is the 
temperature (which is constituted mainly by the internal degrees of freedom of the 
mirror cantilever) and Ax is the separation of two coherent states that are originally 
in a superposition. Note that our experiment is not in the long-time regime where 
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decoherence is characterized simply by a rate. However, the oscillatory term in the 



exponent of Eq. ( 35 i does not affect the order of magnitude and happens to be zero 



after a full period. Assuming that the experiment achieves rj 2 > 1, i.e. a separation 



by the size of a coherent state wavepacket, Ax = y - , the condition that the 
exponent in Eq. p5| l should be at most of order 1 after a full period can be cast in 
the form 

kT 

Q>^=n th , (36) 

where Q = X2m/Jm is the quality factor of the mechanical oscillator. Bearing in 
mind that quality factors of the order of 10 5 — 10 6 have been achieved for silicon 
cantilevers of approximately the right dimensions and frequency, this implies that 
the temperature has to be approximately 3-30 mK. It will be beneficial to perform 
experiments at even lower temperatures to reduce the measurement time, as we will 
explain below. 

Thirdly, the stability requirements for the experiment are very strict. The phase 
of the interferometer has to be stable over the whole measurement time. This means 
that in particular the distance between the large cavity end mirror and the equilib- 
rium position of the small mirror has to be stable to of order A /20N — 0.6 x 10~ 14 m. 

The required measurement time can be determined in the following way. A single 
run of the experiment starts by sending a weak pulse into the interferometer, such 
that on average 0. 1 photons go into either cavity. This probabilistically prepares a 
single-photon state as required to a good approximation. The two-photon contri- 
bution has to be kept low because it causes noise in the interferometer. From Eq. 



(33 1 the width of the revival peak is 2/rjX2M-\/"th- This implies that only a fraction 
~ 1 /ftv^th of the remaining photons will leak out in the time interval of the revival. 
It is therefore important to work at the lowest possible temperature. Temperatures 
below 100 jiK can be achieved with a nuclear demagnetization cryostat. 

Together with the required low value of Qm, the fact that approximately 1% 
of the photons remain after a full period for our assumed loss, and an assumed 
detection efficiency of 70 %, this implies a detection rate of approximately 100 
photons per hour in the revival interval. This means that a measurement time of 
order 30 minutes should give convincing statistics. 

After every single run of the experiment the mirror has to be damped to reset it 
to its initial (thermal) state. This could be done by electric or magnetic fields, e.g. 
following Ref. [48], where a Nickel sphere was attached to the cantilever, whose Q 
could then be changed by 3 orders of magnitude by applying a magnetic field. 

Since the width of the revival peak scales like 1 j\ff , the required measurement 
time can also be decreased by decreasing the temperature below 60 /J.K. Passive 
cooling techniques may be improved. In addition, active and passive optical cooling 
of mirror oscillators has been proposed [49], and implemented experimentally for 
a large mirror (1501] and for small mirrors ||5T1 l52l l53l l54l . Ground state cooling of 
the center of mass motion is achievable and reduces the required measurement time, 
and thus the stability requirements, by a factor of approximately 50. 



Nonclassical States of Light and Mechanics 

4 Entangled states of mechanics and light 



19 



In this last section we turn our attention to non-classical states involving both, light 
and mechanics. In our discussion these will be primarily entangled states, which 
are generated either in steady state under a continuous drive field, or in a regime 
involving short pulses of light. 



4.1 Light mirror entanglement in steady state 

Entanglement of a mechanical oscillator with light has been predicted in a number 
of theoretical studies (SI [55l [58l [59l |60l ET] [62] [63] and would be an intrigu- 

ing demonstration of optomechanics in the quantum regime. These studies, as well 
as similar ones investigating entanglement among several mechanical oscillators 
(^|65]|651|67]|68l|69]|70]|2lJ|72l, explore entanglement in the steady-state regime. 
In this regime the optomechanical system is driven by one or more continuous-wave 
light fields and settles into a stationary state, for which the interplay of optomechan- 
ical coupling, cavity decay, damping of the mechanical oscillator, and thermal noise 
forces may remarkably give rise to persistent entanglement between the intracavity 
field and the mechanical oscillator. 

The simplest example of such a scheme involves an optomechanical system 
driven by one continuous wave laser field. To identify conditions for good optome- 
chanical generation of entanglement we answer a first question that concerns the 
optimal detuning of the driving laser with respect to the cavity field. Given the form 
of the linearized radiation pressure Hamiltonian fig (a + cft)(b + b T ), and the time 
evolution of operators with frequencies A and Q^, we focus on resonant processes 
where A = ±Qm- The first case we analyze is blue -detuning A = —£2m and in 
which we split the interaction in two kind of interactions well known in quantum 
optics: i) beam splitter interaction Jig^b + h.c.) and ii) down-conversion interac- 
tion, hg{ab + h.c). Since the beam splitter term is off-resonant by 2£2y[ and also 
cannnot produce entanglement starting from classical states we drop it and focus on 
the down-conversion term, known to produce bipartite entanglement. Following a 
standard treatment to obtain the covariance matrix in steady state (even analytically 
for this particular case), it can be shown that the logarithmic negativity scales up 
with g as 



<ln 



(37) 



1 +"th 

However, the system is unstable in the "blue-detuned" regime owing to the fast 
transfer of energy from the cavity field to the mirror and an unavoidable bound 
is found g < y^kTm) which in consequence limits E 'jy < ln2. Moreover, any 
thermal quantum n t h > completely destroys the entanglement. One therefore con- 
cludes that the choice of the practical operation regime is dictated by the stability 
of the system. Thus, we move into the "red detuned" which allows for larger g by 
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Fig. 7 (a) Logarithmic negativity E 'jy versus the normalized detuning A / Qm and normalized input 
power P m /P^, (Ff n = 50 mW) at a fixed value of the cavity finesse J? = ,#" = 1.67 x 10 4 ; (b) 
E 'jy versus the normalized finesse J^/J^o and normalized input power Pm/P\ n at a fixed detuning 
A = Qm- Parameter values are X2m/2^ = 10 MHz, J2 = 10 5 , mass m e ff = 10 ng, a cavity of length 
L = 1 mm driven by a laser with wavelength 810 nm, yielding g = 0.95 KHz and a cavity bandwidth 
K = 0.9Qm when = &q. We have assumed a reservoir temperature for the mirror T = 0.4 K, 
corresponding to ~ 833. The sudden drop to zero of E, y corresponds to entering the instability 
region. 



paying the price that, for example at A = Qm the down-converson process is 2Qm 
off-resonant. In this regime analytical results are possible but cumbersome and we 
settled for numerically showing the behavior of E jy in Fig. [7^ as it scales with in- 
creasing input power and varying detunings, and in|7j5 with input power and cavity 
finesse. 

Having concluded that intracavity optomechanical entanglement is attainable, the 
question of detection is to be answered next. As detailed in |57 1 a simple scheme can 
be conceived that consists of a second cavity adjacent to the main one; the second 
cavity output, when weakly driven, does not modify much the first cavity dynamics 
and its output light gives a direct measurement of the mirror dynamics. With homo- 
dyne detection of both cavities and manipulation of the two local oscillators phases 
one can determine all of the entries of the covariance matrix and them numerically 
extract the logarithmic negativity. 

As previously mentioned, in the "red-detuning" regime, the down-conversion 
process is off-resonant and its effect much washed out by the presence of the 
stronger beam-splitter interaction. However, proper detection around the Stokes 
sideband (which carries the photons entangled with the mirror) in the sense de- 
scribed in 2.2 can extract optimal light-mirror entanglement (3). We show this by 
choosing a central frequency of the detected mode Q and its bandwidth T~' and 
computing the bipartite system negativity. The results are shown in Fig. [8^, where 
E, y is plotted versus Q /Qm at five different values of £ = tQm and the other pa- 
rameters similar to the ones used for Fig. [7] If £ < 1, i.e., the bandwidth of the 
detected mode is larger than Qm, the detector does not resolve the motional side- 
bands, and E jy has a value (roughly equal to that of the intracavity case) which 
does not essentially depend upon the central frequency. For smaller bandwidths 
(larger e), the sidebands are resolved by the detection and the role of the central 




frequency becomes important. In particular E jy becomes highly peaked around the 
Stokes sideband £2 = — ^2jvi, showing that the optomechanical entanglement gener- 
ated within the cavity is mostly carried by this lower frequency sideband. What is 
relevant is that the optomechanical entanglement of the output mode is significantly 
larger than its intracavity counterpart and achieves its maximum value at the optimal 
value e ~ 10, i.e., a detection bandwidth T -1 ~ X2m/10. This means that in practice, 
by appropriately filtering the output light, one realizes an effective entanglement 
distillation because the selected output mode is more entangled with the mechanical 
resonator than the intracavity field. 

It is finally important to see what the robustness of the entanglement is with 
increasing temperature of the thermal reservoir. This is shown by Fig. [8J3, where 
the entanglement E jy of the output mode centered at the Stokes sideband is plotted 
versus the temperature of the reservoir at two different values of the bandwidth, the 
optimal one e = 10, and at a larger bandwidth e = 2. We see the expected decay of 
E for increasing temperature, but above all that also this output optomechanical 
entanglement is robust against temperature because it persists even above liquid He 
temperatures, at least in the case of the optimal detection bandwidth e = 10. 



4.2 Entanglement with pulsed light 

An alternative approach to achieving optomechanical entanglement works in a 
pulsed regime. It does not rely on the existence of a stable steady state, such that 
entanglement is not limited by stability requirements. In fact it is possible to op- 
erate in a parameter regime where a stationary state does not exist. This sort of 
optomechanical entanglement can be verified by using a pump-probe sequence of 
light pulses. The quantum state created in this protocol exhibits Einstein-Podolsky- 
Rosen (EPR) type entanglement [6| between the mechanical oscillator and the light 
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pulse. It thus provides the canonical resource for quantum information protocols in- 
volving continuous variable (CV) systems ||73l . Optomechanical EPR entanglement 
can therefore be used for the teleportation of the state of a propagating light pulse 
onto a mechanical oscillator as suggested in I74ll75l . 




Fig. 9 Schematic of the system and the proposed teleportation protocol: (a) A blue detuned light 
pulse (A) is entangled with the mirror (B). (b) A second light pulse (V) is prepared in the input 
state and interferes with A on a beam-splitter. Two homodyne detectors measure P° ut + X v and 
X° ut + P v , yielding outcomes mx and mp respectively. Feedback is applied by displacing the mirror 
state in phase space by a unitary transformation Dx m (mx) Dp m {mp). (c) To verify the success of the 
protocol, the mirror state is coherently transferred to a red detuned laser pulse and a generalised 
quadrature =X 1 ' out cos0 +P|' oul sin 9 is measured. Repeating steps (a)-(c) for the same input 

state but for different phases yields a reconstruction of the mirror's quantum state. 



Let us consider an optomechanical cavity in a Fabry-Prot type setup, cf. Fig. [9]>. 
A light pulse of duration T and carrier frequency C0\ impinges on the cavity and 
interacts with the oscillatory mirror mode via radiation pressure. In a frame rotating 
with the laser frequency, the system is described by the (effective) Hamiltonian HI 

H = Q u Pb + Aa*a+g{b + tf) (a + a 1 ) (38) 

where A — 0) op t — ©i is the detuning of the laser drive with respect to the cavity 
resonance. We assume the pulse to approximately be a flat-top pulse, which has a 
constant amplitude for the largest part, but possesses a smooth head and tail. The 
coupling constant g is then given by 



/ 2KN vh /x 
8 g0 \/(4 2 + K 2 )' 



(39) 
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with A^ph the number of photons in the pulse. It is possible to make either the a 
beam-splitter like interaction (aM + ba^) or the two-mode-squeezing interaction 
(ab + a*b*) resonant by tuning the laser to one of the motional sidebands a) opt ±f2M, 
where the blue (anti-Stokes) sideband (a>\ = <a pt + ^m) enhances down-conversion, 
while the red (Stokes) sideband {coi = 0) op t — &m) enhances the beam-splitter inter- 
action ll76l . In the proposed protocol we make use of both dynamics separately: 
Pulses tuned to the blue sideband are applied to create entanglement, while pulses 
on the red sideband are later used to read out the final mirror state. A similar sep- 
aration of Stokes and anti-Stokes sideband was suggested in ll74l l75ll by selecting 
different angles of reflection of a light pulse scattered from a vibrating mirror in free 
space. 

The full system dynamics, including the dissipative coupling of the mirror and 
the cavity decay, are described by quantum Langevin equations 11771 . which de- 
termine the time evolution of the corresponding operators x m = (b + £> + )/V2, 
p m = —i(b — tf )!^/!. and a,d\ They read 



i m = QMPm, (40) 
p m = -f2 M x m -r M pm-\/2g(fl + fl t ) - \j2r M f, (41) 
— (iA + K)a — iv2gx m — %/2fCa,„, (42) 



a 



where we introduced the (self-adjoint) Brownian stochastic force /, and quantum 
noise a,„ entering the cavity from the electromagnetic environment. Both a,„ and — 
in the high-temperature limit — / are assumed to be Markovian. Their correlation 
functions are thus given by (ai„(t)aj n (t / )) — 8(t — t') (in the optical vacuum state) 
and (f{t)f(t')+f(t')f{t)) = (2n th + l)<5(r-r') (in a thermal state of the mechanics) 

CD- 

We impose the following conditions on the system's parameters. Firstly, we drive 
the cavity on the blue sideband (A = ^2m) and assume to work in the resolved- 
sideband regime (k <C £2m) to enhance the down-conversion dynamics. Note that 
in this regime a stable steady state only exists for very weak optomechanical cou- 
pling [78 1, which poses a fundamental limit to the amount of entanglement that can 
be created in a continuous-wave scheme |5|. In contrast, a pulsed scheme does not 
suffer from these instability issues. In fact, it is easy to check by integrating the full 
dynamics up to time T, that working in this particular regime yields maximal en- 
tanglement, which increases with increasing sideband-resolution Clyi/K. Secondly 
we assume a weak optomechanical coupling g <§; K, such that only first-order inter- 
actions of photons with the mechanics contribute. This minimises pulse distortion 
and simplifies the experimental realization of the protocol. Taken together, the con- 
ditions £2m allow us to invoke the rotating-wave-approximation (RWA), 



which amounts to neglecting the beam-splitter term in 38 Also, we neglect me- 
chanical decoherence effects in this section. We emphasise that this approximation 
is justified as long as the total duration of the protocol is short compared to the 
effective mechanical decoherence time 1 /rj^nth, where Im is the mechanical damp- 
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ing rate and n t h the thermal occupation of the corresponding bath. Corrections to 
this simplified model will be addressed below. 



Based on the assumptions above we can now simplify equations 40 For con- 
venience we go into a frame rotating with Qm by substituting a — > ae ,ilM ', ai„ — > 
ai n e' nM ' and b — > be~' nM '. Note that in this picture the central frequency of a; n is 
located at C0\ — £2m = 0) O pt- In the RWA the Langevin equations then simplify to 

d c = —k& — igb^ — yTKciin, d m = —igd'. (43) 

In the limit g <C K we can use an adiabatic solution for the cavity mode and we 
therefore find 

a(*)« -ilP( t )-^ ain (t), (44) 

b(t) w e Gt b(0) + iV2Ge Gr f de -° s a] n (s) , (45) 



where we defined G = g /k. Equation 45 shows that the mirror motion gets corre 



lated to a light mode of central frequency a>i — Qm (which coincides with the cavity 
resonance frequency ©opt) with an exponentially shaped envelope Um(t) °^ e~ Gt . Us- 
ing the standard cavity input-output relations a ou , = a ln + \/2Ka allows us to define 
a set of normalised temporal light-modes 



,-2Gt 



de - Q, a in (f ) , A out = J ^ _ - J de Gt a out (t ) , (46) 



which obey the canonical commutation relations [A,-,Aj ] = 1. Together with the def- 
initions Z?; n = b(0) and B out = b{x) we arrive at the following expressions, which 
relate the mechanical and optical mode at the end of the pulse t = T 



Aout = -e Gx A m - We^-lBl , B out = e Gx B m + iy/ ' - ■ IaL . (47) 



By expressing equations 47 in terms of quadratures X' m — (B,- + Bl)/v2 and X[ = 
(A; +Aj)/v / 2, where i e {in, out}, and their corresponding conjugate variables, we 
can calculate the so-called EPR-variance Aepr of the state after the interaction. For 
light initially in vacuum (AX™) 2 — (AP™) 2 = \ and the mirror in a thermal state 
(AX™) 2 = (AP™) 2 = n + \, the state is entangled iff CD 

A EPR = [A (XT + Pr 1 )} 2 + [A (PT + XT)} 2 (48) 
= 2(« + l)(V-v / e 2 ''-l) 2 <2 ! (49) 

where r = Gt is the squeezing parameter and hq the initial occupation number of the 
mechanical oscillator. Note that in the limit of large squeezing r>lwe find that the 
variance Aepr ~ (no + l)e~ 2r /2 is suppressed exponentially, which shows that the 
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created state asymptotically approximates an EPR-state. Therefore, this state can be 



readily used to conduct optomechanical teleportation. Rearranging 48 we find that 
the state is entangled as long as 

r > l taB0) (50 ) 



2 V 4 (»o + l)/ 2 

where the last step holds for no — > °° ■ This illustrates that in our scheme the require- 
ment on the strength of the effective optomechanical interaction, as quantified by 

2 

the parameter r = scales logarithmically with the initial occupation number no 
of the mechanical oscillator. This tremendously eases the protocol's experimental 
realization, as neither g nor T can be arbitrarily increased. Note that no need not be 
equal to the mean bath occupation n t h, but may be decreased by laser pre-cooling to 
improve the protocol's performance. 

To verify the successful creation of entanglement a red detuned laser pulse (A = 
£2m) is sent to the cavity where it resonantly drives the beam-splitter interaction, 
and hence generates a state-swap between the mechanical and the optical mode. It 
is straightforward to show that choosing A = £2m leads to a different set of Langevin 



equations which can be obtained from 43 by dropping the Hermitian conjugation (f ) 
on the right-hand-side. By defining modified mode functions 0£^ ollt j = 0£ O ut(in) m ^ 
corresponding light modes A' in ^ out ^ one obtains input/output expressions in analogy 
to ((47) 



Km = -e~ Gx A' m + iVl-e- 2G *B in , B out = e - Gx B m - iy/l-e'^A^. (51) 

The pulsed state-swapping operation therefore also features an exponential scaling 
with Gt. For Gt — > °° the expressions above reduce to A' ou{ = — /B; n and B out = iA' m , 
which shows that in this case the mechanical state — apart from a phase shift — is 
perfectly transferred to the optical mode. In the Schrdinger-picture this amounts to 
the transformation |<p) m | y)\ — > | V)m|<p)i. where <p and yf constitute the initial state 
of the mechanics and the light pulse respectively. The state-swap operation thus 
allows us to access mechanical quadratures by measuring quadratures of the light 
and therefore to reconstruct the state of the bipartite system via optical homodyne 
tomography. 

As we have shown above, pulsed operation allows us to create EPR-type entan- 
glement, which forms the central entanglement resource of many quantum informa- 
tion processing protocols [73 1. An immediate extension of the proposed scheme is 
an optomechanical continuous variables quantum teleportation protocol. The main 
idea of quantum state teleportation in this context is to transfer an arbitrary quantum 
state | y/jn) of a travelling wave light pulse onto the mechanical resonator, without 
any direct interaction between the two systems, but by making use of optomechani- 
cal entanglement. The scheme works in full analogy to the CV teleportation proto- 
col for photons l80l IHTl and, due to its pulsed nature, closely resembles the scheme 
used in atomic ensembles [82, 83 1: A light pulse (A) is sent to the optomechanical 
cavity and is entangled with its mechanical mode (B) via the dynamics described 
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above. Meanwhile a second pulse (V) is prepared in the state |y/m), which is to 
be teleported. This pulse then interferes with A on a beam-splitter. In the output 
ports of the beam-splitter, two homodyne detectors measure two joint quadratures 
P° ut +X v andX° ut +P V , yielding outcomes mx and mp respectively. This constitutes 
the analogue to the Bell-measurement in the case of qubit teleportation and effec- 
tively projects previously unrelated systems A and V onto an EPR-state [84|. Note 
that both the second pulse and the local oscillator for the homodyne measurements 
must be mode-matched to A after the interaction, i.e., they must possess the iden- 
tical carrier frequency as well as the same exponential envelope. The protocol is 
concluded by displacing the mirror in position and momentum by mx and mp ac- 
cording to the outcome of the Bell-measurement. This can be achieved by means of 
short light-pulses, applying the methods described in ll20l [85l . After the feedback 
the mirror is then described by ll73l 

X£ = XT + Pr +X V =X V + (V - y/e^l) PC - ff) , (52) 

Pm = PT +XT +Pv=Pv+ (V - V^l) (C - X( a ), (53) 

which shows that its final state corresponds to the input state plus quantum noise 
contributions. It is obvious from these expressions that the total noise added to both 
quadratures (second term in [52] and [53] respectively) is equal to the EPR-variance. 
Again, for large squeezing r> 1 the noise terms are exponentially suppressed and 
in the limit r4«i, where the resource state approaches the EPR-state, we obtain 
perfect teleportation fidelity, i.e., X^ 1 — X v and P^ n — P v . In particular this operator 
identity means, that all moments of X v , P v with respect to the input state | Yin) Wl H be 
transferred to the mechanical oscillator, and hence its final state will be identically 
given by 

We found that in the ideal scenario the amount of entanglement essentially de- 
pends only on the coupling strength (or equivalently on the input laser power) and 
the duration of the laser pulse and that it shows an encouraging scaling, growing 
exponentially with Gt. This in turn means that the minimal amount of squeezing 
needed to generate entanglement only grows logarithmically with the initial me- 
chanical occupation no. In a more realistic scenario one has to include thermal noise 
effects and effects of counter-rotating terms. Including the above-mentioned per- 
turbations results in a final state which deviates from an EPR-entangled state. To 
minimise the extent of these deviations, the system parameters must obey the fol- 
lowing conditions: 

1 . K <C X2m results in a sharply peaked cavity response and implies that the down- 
conversion dynamics is heavily enhanced with respect to the suppressed beam- 
splitter interaction. 

2. g < K inhibits multiple interactions of a single photon with the mechanical mode 
before it leaves the cavity. This suppresses spurious correlations to the intracavity 
field. It also minimises pulse distortion and simplifies the protocol with regard to 
mode matching and detection. 
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3. gx ^> 1 is needed in order to create sufficiently strong entanglement. This is due 
to the fact that the squeezing parameter r = (g/K)gx should be large, while g/x 
needs to be small. 

4. n t hlMT <C 1, where is the thermal occupation of the mechanical bath, assures 
coherent dynamics over the full duration of the protocol, which is an essential 
requirement for observing quantum effects. As the thermal occupation of the 
mechanical bath may be considerably large even at cryogenic temperatures, this 
poses (for fixed Im and n t h) a very strict upper limit to the pulse duration T. 

Note however that not all of these inequalities have to be fulfilled equally strictly, 
but there rather exists an optimum which arises from balancing all contributions. It 
turns out that fulfilling (4) is critical for successful teleportation, whereas (l)-(3) 
only need to be weakly satisfied. Taking the above considerations into account, we 
find a sequence of parameter inequalities 

«thJM< - <g< K«f2 M , (54) 

T 

which defines the optimal parameter regime. Dividing this equation by Tm and tak- 
ing a look at the outermost condition n t h <C Q, where Q — X2m/Fm is the mechanical 
quality factor, we see that the ratio Q / defines the range which all the other pa- 
rameters have to fit into. It is intuitively clear, that a high quality factor and a low 
bath occupation number, and consequently low effective mechanical decoherence, 
are favourable for the success of the protocol. Equivalently, we can rewrite the occu- 
pation number as «,h = ksT^h/hQu and therefore find k^T^^/h <C Q ■ £2m, where 
now the Q ■ /-product (/ = £2m/2tz) has to be compared to the thermal frequency 
of the bath. Let us consider a numerical example: For a temperature 7b a th ~ lOOmK 
the left-hand-side gives &B7bath/ft ~ 2ft- 10 9 Hz. The Q- /-product consequently has 
to be several orders of magnitude larger to successfully create entanglement. As 
current optomechanical systems feature a Q ■ /-product of 2n ■ 10 11 Hz and above 
[86, 87, 88, 89 1, this requirement seems feasible to meet. 



5 Conclusion and Outlook 

The selection of protocols presented in the present chapter clearly demonstrate the 
feasibility — and the stringent requirements — for quantum state engineering in 
mesoscopic mechanical systems. At the time of writing this book chapter state of 
the art experiments achieve ground state cooling, but manifest quantum effects in the 
dynamics of optomechanical systems have yet to be demonstrated. From the discus- 
sion given above it should be clear that a necessary condition for observing quantum 
effects is a sufficiently large product of mechanical quality factor and frequency, 
Q£2m > kgT /h. What exactly "large" means in this context depends crucially on 
the protocol to be implemented. The goal of making the quantum regime accessible 
for mechanical systems thus has to be approached from both sides: Experimentally, 
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by developing optomechanical systems with a sufficiently large QQm -product; and 
theoretically, by developing schemes which are not too demanding regarding the 
magnitude of this number. Once those two ends meet this will mark the birth of a 
new field of research, quantum optomechanics. 
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